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Specific lieat and low-lying excitations in the mixed state for a type II superconductor 
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Low temperature behavior of the electronic specific heat C{T) in the mixed state is by the self- 
consistent calculation of the Eilenberger theory. In addition to 7r-term (7 is a Sommerfeld coeffi- 
cient), C{T) has significant contribution of T^-term intrinsic in the vortex state. We identify the 
origin of the T^-term as (i) V-shape density of states in the vortex state and (ii) Kramer-Pesch 
effect of vortex core shrinking upon lowering T. These results both for full-gap and line node cases 
reveal that the vortex core is a richer electronic structure beyond the normal core picture. 

PACS numbers: 74.25.0p, 74.25. Bt, 74.25.Jb 



There has been considerable attention focused on the 
properties of a type II superconductor under an applied 
field, ranging from conventional to exotic pairing super- 
conductors, such as cuprates. Interest arises not only 
from a fundamental physics point of view, but also from 
an application point of view, such as the quest to achieve 
a higher critical current, etc. Since practical applica- 
tions of superconducting technology still predominantly 
use "conventional" compounds, such as A- 15 compounds; 
NbsSn or NbsAl as superconducting magnets, it is of con- 
siderable importance to thoroughly understand the fun- 
damental properties of the mixed state in an ordinary s- 
wave pairing state. In spite of a long historjti, it is only 
recently that we are discovering the detailed electronic 
structures of the mixed state experimentalljiSii^ and 
theoretically^. The latter has become possible because 
of combined efforts using both analytical and numerical 
computations based on a microscopic model. 

The specific heat can be expanded as 
C = 7T -I- C2T^ -!-••• at a low temperature T in 
the mixed state. The first term jT is the electronic 
contribution proportional to the zero-energy density of 
states (DOS) N{E = 0). The 7-term has been well 
studied for the purpose of determining the electronic 
state. For example, we can investigate the gap topology 
using the magnetic field B dependence; 7 ex i? for 
s-wave full gap superconductor and 7 ex \/B for d-wave 
superconductor with line nodes in the superconducting 
gap^-^. However, so far, the contribution of the T^-term 
to C has not been discussed^. The purpose of this letter 
is to quantitatively estimate this T^-term by microscopic 
calculation, and to clarify the existence and origin of 
the r^-term. For this purpose, we also investigate the 
spectrum of the DOS N{E) in the mixed state and 
T-dependence of the vortex structure. From this analysis 
we find that the T^-term in C comes from the electronic 
features intrinsic in the vortex state, i.e., V-shape DOS 
spectrum as seen later and the shrinkage of the core 
radius with decreasing T, namely, the Kramer-Pesch 
(KP) effecliSiifliii. When the vortex core radius keeps 
shrinking, N{E = 0) is expected to acquire an extra 
T-dependence. The KP effect is one of the consequences 



of the fact that the vortex can not be pictured as a 
simple rigid core filled by normal electrons. The physics 
of vortices at low T is much more interesting than the 
"normal core model" . 

Recently, several careful experiments on PbMogSgi^, 
Nb77Zr23iS, /3-pyrochlor CsOszOgi^ and RbOsjOgi^ in- 
dicate a hint of a substantial deviation from the simple 
7T behaviors even at very low T. These materials are 
regarded as "conventional" superconductors. 

Here we are going to numerically calculate C{T), N{E) 
and the vortex core size p by solving the quasi-classical 
Eilenberger equation self-consistently with greater nu- 
merical precision than has hitherto been reported. Do- 
ing this is crucial for examining the distinct, but subtle 
low T behaviors of these quantities and allows us to find 
the internal relationship between them. If lower numer- 
ical precisions are employed computational results yield 
only rounded features that lead to the normal rigid core 
model. We assume an isotropic s-wave superconductor 
in the clean limit. To our knowledge, there has been 
no previous study on C{T) and N{E) to this level of 
accuracjiii2ii4. 

The quasi-classical theory is valid for the case when 
/cfC ^ Ij which is satisfied for almost all type II super- 
conductors, fcp is the Fermi- wave number and ^ is the 
BCS coherence length which is our units of the length 
scale. We introduce the pair potential A(r), the vector 
potential A(r) and the quasi-classical Green's functions 
g{iujn,'r,0), f{ioJn,r,9,) and f''{iujn,r,9), where r is the 
center of mass coordinate of the Cooper pair. The di- 
rection of the momentum, k = k/ |k|= (cos0,sin6'), is 
represented by the polar angle 6 relative to x-direction. 
The Eilenberger equation is given by 
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where g = [1 — f'^f]'^, Rc^ > 0, vp = wpk is the Fermi 
velocity, and (j)o is a flux quantumSi. The applied field 
H is along the z-direction. In the symmetric gauge, the 
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FIG. 2: (Color online) Energy dependence of the total density 
of states N{E) is calculated at T/Tc = 0.01. Plotted data are 
for B/Bo = 0, 0.3, 0.4, 0.5, 0.6 from bottom to top at S = 0, 
showing characteristic V-shape cusp structures near the Fermi 
level E = Q. 



FIG. 1: (Color online) (a) Free energy difference Fg. — Fn 
and (b) specific heat C{T)/T as a function of T/Tc for var- 
ious fields B/Bq = 0, 0.3, 0.4, 0.5, 0.6 from bottom to top. 



vector potential is written as A(r) = (1/2)H x r + a(r), 
and the internal field h(r) is given by h(r) = V x a(r). 

Equation Q) can be self-consistently solved by numer- 
ical calculation, by considering the self-consistent condi- 
tions for A(r) and a(r); 



at an energy E are given by 



A{r) = NoVo2nT J2 ^/(ic^n,r, 0'), (2) 

J(r)=-^2.r5:y -.5(-„,r,.), (3) 

where j(r) = V x V x a(r) and A'o is the density of states 
at the Fermi level in the norm al state. T he cut-off en- 
ergy is set as wc = 20re. n = \/7C(3)/72(Ao/rc)«;GL. C 
is Riemann's zeta function. The Ginzburg-Landau pa- 
rameter is chosen as kgl = 9.0. Aq is the uniform gap 
at T = 0. 

Using the self-consistent solutions g, /, /^ A and a, 
the free energy difference is given by 
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where Fs{Fn) is the free energy for the superconducting 

(normal) state, and (• • •)r indicates the spatial average. 

From Fs, the specific heat coefficient C/T is obtained by 

C/T = -d^Fs/dT^ . The local density of states (LDOS) 
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where g is calculated by Eq. (Q) with ia;„ -^ E + iri. We 
set r] = O.OIAq. The total DOS N{E) is the spatial av- 
erage of the LDOS, i.e. N{E) = {N{E,r))^. The self- 
consistent calculation is performed within the vortex lat- 
tice unit cell, which is divided into 81 x 81 mesh points. 
We assume that the vortices form a triangular lattice. 

The free energy difference Eg — F^ and the specific heat 
C{T)/T at low T are shown in Fig. ^for various fields 
B, which are normalized by Bq = (/)o/C^- We confirm 
that C{T)/T in our numerical calculation at i? = is 
identical with the BCS result. It is clearly seen that for 
B ^0, C{T)/T exhibits a T-linear behavior at lower T 
{T/Tc < 0.2) where C{T) for B = vanishes exponen- 
tially. Thus the specific heat in the mixed state can be 
expressed as 



C{T) _^ , T_ 

^^-7 +«c^ 




(6) 



at low T, where 7' = 7/7,1 = N{E = 0)/Nq 
and 7n = 2tt'^Nq/3. We also confirm that 

7'(_B) = \imT^QC/{jnT) shows a _B-linear dependence 
expected for the s-wave pairing state^. It is remarkable 
to notice that C/T is not constant at finite T, indicating 
a significant contribution from the acr-tcrm. When we 
estimate 7' from C/T, the experimental data has to be 
carefully extrapolated toward T — ^- 0, in order to remove 
the acT-term contribution. 

The contribution of ac^^-term in Eq. © also reveals 
important information on the electronic states. In order 
to analyze the origin of the acT-term, we calculate the 
DOS spectrum and estimate the KP effect. The total 
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FIG. 3: (Color online) Spectral evolution A'^(i5,rNN), where 
tnn is along the nearest neighbor vortex direction of the 
hexagonal vortex lattice. The vortex center is situated at 
ItnnI = 0. A''(0, 0) is truncated over 20. 



DOS N{E) is shown in Fig. Elat our lowest r(= O.OlTc). 
We clearly see the V-shape DOS with a cusp structure 
for B ^ 0. The low energy part of N{E) can be fitted as 



N{E) _ N{E^O,T) 
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As we show later more quantitatively, N(E — 0, T) and 
a_E correspond to T-linear and acT"^ contributions, re- 
spectively, to the specific heat C{T). From Fig. [3 we 
also see that N{E = 0) ^^ j increases with increasing B 
and that as decreases as B increases, so that the linear 
behavior of \E\ is unchanged. 

The i?-dependence of N{E) is related to the spec- 
tral evolution of the LDOS N{E, r) around the vortex, 
which is shown in Fig. 13 It is clear that the quasi- 
particle bound states have a characteristic dispersion re- 
lation \E\ c>c r as a function of the radius r = |r| from 
the vortex center, and form ridges extending outwards 
from the core. The zero-energy peak at the vortex cen- 
ter gives rise to N{E = 0), and the finite energy low- 
lying bound states near the core ultimately gives the 
|ii^|-dependent term in N{E). A simple understanding 
of Eq. (Q at lower B or lower E is as follows: The 
total DOS can be evaluated by integrating LDOS spa- 
tially. Assuming circular symmetry around the core, 
N{E) = N{E ^ 0) + 2TrJ^S{\E\ - /3r)rdr per unit 
length along the field direction. The first term N{E = Q) 
mainly comes from LDOS N{E — 0, r) at the vortex 
core r = (see the peak at i? = in Fig. |21l. The sec- 
ond term comes from the quasi-particle spectral weight 
whose spatial trajectory \E\ = /3|r|, is indicated by the 
ridges in Fig. O This gives rise to V-shape total DOS: 
N{E) = N{E == 0) + 2TT\E\/f3^. These rich features are 
distinct from the simple normal core model, which yields 
a featureless dispersion around the core. Some of the 
features of LDOS are directly observed in NbSe^^, VaSi^ 
and YNiaBaGi via STM-STS experiments. 

We now estimate quantitatively the contribution from 
the KP effect to C{T)/T. Here we define the core radius 
p from the slope of A(r), i.e., p = Ao(9A(r)/9r)~^ at the 
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FIG. 4: (Color online) (a) Temperature dependence of the 
vortex-core radius p and (b) Temperature dependence of the 
zero-energy density of states N{E — 0, T) are plotted for 
B/Bo = 0, 0.3, 0.4, 0.5, 0.6 from bottom to top. In Fig. (a) 
data for B — and normal state are not shown. Lines are 
guides for the eye. 
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FIG. 5: (Color online) Specific heat C{T)/T for B/Bo = 0.3 
(solid circles) with the estimated contributions LdasT/Tc, 
aKT/Tc, and (1.9aB + aK)T/n. 



vortex core along the nearest neighbor vortex direction. 
As shown in Fig. Efa), p decreases almost linearly on 
lowering T for various values of B, which is a confirma- 
tion of the KP effect. It is noted that the limit of p when 
T ^ is finite in the vortex lattice case, while p — > at 
T — + in the single vortex casein. Since the zero-energy 
DOS due to the bound state around the vortex core is 
related to the core radius p, N{E = 0,T) decreases as T 
decreases, as shown in Fig. EJb). This is approximately 
fitted by 
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Combining Eqs. 

we can express 

N{E,T)/No = Y + aKT/T, 
relation C{T)/T = {2/T) J° 
using the Fermi 

/(^,T) and Ao/T, 

C(r)/(7„r) = 7' -I- (aK + 1.9aE)T/T,. 
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0, 



and 

the DOS as 

'raE\E\ / Ao . Through the 

' dEEN{E, T)df{E, T)/dT 

distribution function 

= 1.76, we obtain 

It is now 



clear that the aJT term in C{T)/T [Eq. ®] is affected 
by the KP effect term ax and i?-linear coefhcient aE 
equally. To examine this relation, we analyze the data 
at B = 0.3 as an example. In Fig. |31we display C{T)/T 
shown in Fig. [Hb), together with 7' + LOafiT/Tc 
and 7' + (aK + 1.9a£;)T/Tc. It can be seen that these 
two contributions uk and as account for almost all 
C{T)/T behavior at low T. That is, the KP term uk 
and i?-linear term as are both equally important for 
understanding the total specific heat quantitatively. 

We also study the d-wave pairing case with a line node, 
where A(r) -^ A(r)\/2cos26' in Eq. ^^'"^ , and obtain 
the same results as for the s-wave case: V-shape DOS, 
KP effect and a^r-term in C/T^^. Note that C{T)/T 
behaviors in Fig. ^b) arc qualitatively identical to the 
data of Sr2Ru04 with a line node gap^^. 

Our results of OcT-term in C /T and V-shape DOS can 
be obtained after performing an exact calculation using 
the self-consistently obtained pair potential A(r). When 
C{T)/T and N{E) are calculated by solving the same 
Eilenberger equation within the Pesch approximation, it 
is found that there is neither a linear Q!£;|-E| term in the 
DOS nor a linear term acT in the specific heat (see Figs. 
1(b) and 1(d) inRef. 14). 

We discuss briefly the relationship between the V- 
shape DOS and physical quantities. The |_E|-linear func- 
tional dependence in the V-shape DOS N{E) happens to 
be the same as the DOS N{E) ex \E\ for the line node 
gap structure at a zero field. From the |_E|-linear de- 
pendence, we can evaluate the power law T-dependence 
for various physical quantities by simple power-counting, 
such as specific heat {C{T) ex T^), nuclear relaxation 
time (Ti(T)-i oc T^), thermal conductivity (k(T) cx T^), 
and ultrasound attenuation. These behaviors are used to 
identify the line-node of the superconducting gap. In the 
mixed states, these power law components appear due to 
the V-shape DOS both for the full-gapped and line-node 
cases, in addition to the contributions by the zero-energy 
DOS N{E = 0). Therefore, in the experiment when a 
magnetic field is applied, we can not simply assign the 
origin of the power law behavior as a line node, since the 
origin may be the V-shape DOS due to the vortex states. 
The Ti-behavior Ti{T)~^ ex T^ under magnetic fields is 



a consequence of the spatial average. If we observe Ti 
outside of the vortex core by using a site selective NMR 
techniqueii, we can unambiguously identify the signal 
due to the line node. 

In summary, we have demonstrated that the elec- 
tronic specific heat at low T obeys the general law 
C{T)/T = Y + ttcT/Tc in the mixed state of full-gapped 
superconductors. The contribution of second term 
UcT/Tc, caused by the electronic features intrinsic in 
the vortex state, is significant in C /T and observ- 
able. Correspondingly we have found the V-shaped DOS 
N{E) = N{E ^ 0)-FaB|£;|/Ao in the vortex states. The 
KP effect, shrinking vortex core radius with decreasing T, 
contributes to the T-dependence of the zero-energy DOS; 
N{E = 0,r)/A^o = 7' + dKT/T^. We have shown that 
ftc ^ 1.9aB -I- ax, which shows the contributions from 
both the V-shape DOS and the KP effect are responsi- 
ble for the QfcT/Tc term in C{T)/T. We have also cal- 
culated for extreme anisotropic-gap cases, namely nodal 
gap cases, and found that the V-shape DOS and KP ef- 
fect are not altered -'^^ . We also pointed out some cautions 
when attempting to identify the gap topology experimen- 
tally, since the line node gap yields the same V-shape 
DOS even in the full gap under a field. 

We expect that careful measurement of C{T) for clean 
type II superconductors should confirm our prediction. 
The V-shape DOS is directly observable by tunneling ex- 
periments. These characteristics are found to be deeply 
rooted in the low-lying excitations around a vortex core, 
which goes beyond the rigid normal core picture, reveal- 
ing a richer electronic structure in the vortex core. We 
emphasize that the sample quality is important; Impurity 
effects nulhfy the T-linearity C{T)/T. 
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